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Synopsis 

We study both experimentally and theoretically the rheological behavior of 
isotropic bidisperse suspensions of noncolloidal particles in yield stress fluids. 
We focus on materials in which noncolloidal particles interact with the suspend- 
ing fluid only through hydrodynamical interactions. We observe that both the elastic 
modulus and yield stress of bidisperse suspensions are lower than those of monodisperse sus- 
pensions of same solid volume fraction. Moreover, we show that the dimensionless yield 
stress of such suspensions is linked to their dime nsionless elastic modu lus and to their solid 
volume fraction through the simple equation of 



Chateau et al 



(|2008[ ). We also show that 



the effect of the particle size heterogeneity can be described by means of a packing model 
developed to estimate random loose packing of assemblies of dry particles. All these obser- 
vations finally allow us to propose simple closed form estimates for both the elastic modulus 
and the yield stress of bidisperse suspensions: while the elastic modulus is a function of the 
reduced volume fraction (j)/4>m only, where 0m is the estimated random loose packing, the 
yield stress is a function of both the volume fraction (jj and the reduced volume fraction. 

I Introduction 

Many industrial processes such as concrete casting, drilling muds, food-stuff transport ... and 
natural phenomena, such as slurries, lava flows ... involve suspensions of polydisperse particles 
suspended in a non-Newtonian fluid. Knowing and predicting the rheological properties of such 
suspensions is thus a major issue of both industrial materials mix design and science of defor- 
mation and flow of materials. For instance, it is well known that the grading of aggregate is 
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one of the main factor influencing the hardened concrete stren gth; on the ot her hand, it has 



Nevillel (|l98lh ]. Consequently, 



a large impact on the workability of concrete in its fresh state 
a lot of work has been devoted to predict the influence of the particle size distribution on the 
concrete overall prop erties, both in the hardened and fresh state, and then to elaborate rational 



mix design methods |de LarrardI (Il999l )]. As the rheological properties of a suspension (viscosity, 
yield stress, ...) are increasing functions of the solid particle volume fraction, it is necessary to 
lower this effect in order to design concrete mixtures containing the maximum amount of solid 
particles possible that can be transported, placed and finished easily. For this purpose, it is 
better to use aggregate with a distributed grading rather than particles of similar sizes. 

The behavior of polydisperse, in particular bidisperse, Newtonian suspensions has received 
large interest from rheologists and several works, either theoretical or experimental, have been 
devoted to this subject in the last few decades. Different theoretical relationships aiming at pre- 



dicting the viscosity of polydisperse suspensions as a function of the solid vo ^ 
of the particle size di s tribut i on have been proposed in the literature 



ami); 



Storms et al. 



FarrisI (jl968|) 



ume fraction 6 and 



Chong et al. 



Gondret and PetitI (jl997l )]. In the three last works, the influence 



of the particle size distribution is taken into account through the concept of maximum packing 
fraction even if this concept is not always rigorously defined. It seems quite natural to define the 
maximum packing fraction of a suspension as the volume fraction for which the rheological prop- 
erties tend to diverge. Nevertheless, the determination of the value of the sohd volume 
fraction at v^hich the rheological properties diverge, and its link to packing models, 
pose problem. In th e case of monodisperse spheres, it seems to be equal to 0.5 7 for isotropic 



disper sions of particles 



Mahaut et al. 



(l2008al )]. to 0.605 for anisotropic dispersions 



Ovarlez et al. 



(120061)1 w h ile th e random close packing of the dry particles is equal to 0.64. 



FarrisI (j 19681 ) has investigated the relative viscosity of bimodal suspensions in the framework 
of a rigorous theoretical approach, valid only if the coarse particle size is much larger than the 
fine particle one (i.e. when A, the coarse to fine particle ratio, is greater than 10). When this 
condition is fulfilled, one can consider that the coarse particles interact with the Newtonian fine 
particle suspension, and thence, the relative viscosity of the bimodal suspension is equal to the 
product of the relative viscosities of each monodisperse suspension. Of course, Farris approach 
can be generalized to polydisperse suspensions or to non-Newtonian suspei iding fluids provided 



that estimates of the monodisperse suspension's overall properties exist. 



Stovall et al 



(119871 ) 



generalized Farris approach by taking into account interactions between the particles of different 
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size using the random close packing's model of iStovall et al. 



(j 19861 ) in order to address prob- 



lems where scale separation between particles of different size is not possible. Nevertheless, it is 
believed that this work lacks of rigor because Farris reaso ning is used by the autho rs in a situa- 



Phan-Thien et al. 



(| 19971 ) developed a 



tion where particle size separation is not possible. Later, 
multiphase model for polydisperse suspensions. The idea consists in starting from the homoge- 
neous linear suspending material (a Newtonian fluid or a Hookean solid) and in introducing the 
particles by infinitesimal volume fraction in the framework of an iterative process (differential 
scheme). The first step actually corresponds to the Einstein approach to the overall properties 
of a dilute suspension. Next steps consist in removing a small volume fraction of the suspending 
overall medium and in replacing it by the same volume of particles. Of course, such a process is 
rigorous only if particle size separation is possible at each step. Then, this model is only valid for 
the l ow solid volume fr a ctions and hence is not applicable to concentrated suspensions. By the 



way, 



Phan-Thien et al. 



(|l997l ) found that the values of the rheological properties estimated by 
this model are much smaller than those experimentally measured on polydisperse concentrated 
suspensions. 



Bes ides theoretical works, exper imental studies hav e also 



(11992) 



Probstein et al 



(119941 ) and 



seen performed. 



Shapiro and Probstein 



Chang and Powelll (jl994l ) used different types of viscometers 



to measure the steady viscosity of bidisperse suspensions. The same general trends are observed 
in these works. For a given value of the solid volume fraction, the viscosity of a bidisperse sus- 
pension is lower than the viscosity of a monodisperse suspension with same solid volume and 
the viscosity is a decrea sing function of the maxi mum packing fraction of the particle size dis- 



tribution. Interestingly, IChang and Powelll (1 19941 ) also showed that the dimensionless viscosity 
rj{(p)/rj{0) of such suspensions is basically a function of (p/cjym only, where is the particle vol- 



ume fraction and (J )m is claimed to be the maximum packing fraction of the particle mixture. 



Gondret and Petit 



(|l997l ) have measured the finite frequency viscosity of bidisperse suspensions. 
They also observed that the particle size distribution influences the dynamic viscosity of bidis- 
perse suspensions, and that the lower maximum packing density distributions correspond to the 
lower dynamic viscosities. 

>From these studies, it clearly appears that the viscosity of a polydisperse concentrated 
suspension is closely related to the maximum packing fraction of the suspended particles. Ac- 
cordingly, the vi scosity of a Newtoni an suspension can be controlled by optimizing the particle 



size distribution 



Servais et al. 



(12002)]. 
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All these studies focused on the influence of the polydispersity on the viscosity of Newtonian 
suspensions while polydisperse non-Newtonian suspensions have been poorly studied. To our 
knowledge, this problem has been studied only by few a uthors and the i r studie s pro vide extremely 



dispers ed results. However, as explained in detail by 



Mahaut et al 



(2008a) and 



Mahaut et al. 



(|2008bl ). the data reported in these works do not correspond to homogeneous suspensions of 



particles interacting only mechanically with the suspending fluid, whic h is the subject we are 



Geiker et al. 



(12003) 



interested in, and then are not applicable to a wide range of materials, 
studied experimentally the effect of coarse particles volume fraction on the rheologi- 
cal properties of self compacting concretes; they studied suspensions of polydisperse 
particles of various shapes but of same grading. They assumed that the effect of 
aggregates on fresh concrete rheological properties can be studied by looking to 
concrete as a suspension of coarse particles (the aggregates) in a yield stress fluid 
(the mortar): coarse particles interact only by non-Newtonian hydrodynamic forces. 
Both the fresh concrete a nd the mortar behaviors are described by a Bingham law^. 
While 



Geiker et al. 



(|2002l ) find that the dimensionless rheological properties for the 



different studied particles have very different shapes v^hen plotted vs. the reduced 
volume fraction (p/cpm^ close examination of the results show that divergence actually 
occurs at different (p/cpm < 1 for particles of different shape, in contradiction with 
the usual definition of the maximum packing fraction in rheology. In this study, 
was indeed defined as the random close packing of the dry granular assembly, which 
may not be relevant for comparing the rheological data. Moreover, the concrete 
was tested in a large gap coaxial cylinder rheometer and its rheological properties 
were estimated from the steady-state flow. It is well known that shear ind uced mi- 



gration of particles is likely to occur in this situation [Ovarlez et al 



(|2006h ]. Then, 



the tested materials may not be homogeneous w hen the syste n i is a t steady-state 



Geiker et al. 



(120021 ) can be used 



and it is not clear that the experimental data of 
to estimate the rheological properties of a well defined material (i.e. a material ho 
mogeneous at the coarse particles scale with a well defined distributio n of particles 



within the tested volume). The fact that the yield stress estimates of 



Geiker et al. 



((20021) are much larger than estimates reported by other authors or in this work 



als o casts doubt abou t the v alidity of their procedure 



Ancey and JorrotI (120011 ) have experimentally studied the influence of adding 
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noncolloidal and non-Brownian particles within a clay dispersion on the value of the 
yield stress of the suspension. The ratio of the particle size was large enough so that 
the clay dispersion can be considered as homogeneous at the coarse particles scale. 
The authors have chosen to measure the yield stress of the suspension by means 
of a slump test in order to avoid the classical problems encountered with Couette 
rheometer (migrat i on of p articles, locali z ation of the shear rate, anisotropy of the 



material) 



CoussotI (120051) 



Ovarlez et al. 



(j2006l )]. They showed that for well-graded 



materials, the monodisperse suspension yield stress does not depend on the particle 
characteristics (diameter, material) and that the relative yield stress diverges when 
the solid volume fraction value tends toward that of the maximum packing density. 
When the coarse particles are polydisperse, the value of the maximum packing 
density depends of the size distribution of the particles. Similarly to what is observed 
for Newtonian suspensions, the yield stress diverges for values of the solid volume 
fraction depending on the particle size distribution. The dimensionless yield stress 
Tc{'P)/tc{0) measured for various mixtures of bidisperse particles then seem to collapse 
when plotted vs. the reduced volume fraction <j)/4>m- They also observed that for 
low reduced solid volume fraction, the yield stress can be a decreasing function of 
the solid volume fraction of the coarse particle. This effect was ascribed by the 
authors to a depletion phenomena. In the closeness of coarse particles, the clay 
particles are expelled from the suspending fluid. This efl'ect induces an increase of 
the clay particle concentration far from the coarse particles which are embedded in 
a shell of viscous material (pure water) where the yield stress is naught. Then, the 
coarse particles can not contribute to the overall yield stress and behave as empty 
pores. Obviously, the observed decrease of the suspension yield stress originates 
from physicochemical effects. We recall that such efl"ects are beyond the scope of 
this paper. 

In this study, we are interested in situations where bidisperse mixtures of noncolloidal par- 
ticles are dispersed in a yield stress fluid. It is recalled that yield stress fluids have a solid 
viscoelastic behavior below a yield stress; above this yield stres s, they behav e as liquids and 



their flow behavior is often well fltted to a Herschel-Bulkley law 



LarsonI (jl999l )]. The influence 



of monodisperse particle s on t 



by 



Mahaut et al. 



pnnSal ) and 



le behavior of yield st ress fluid has been ad dressed experiinental 



Mahaut et al\ (|2008bl ) , and theoretically by 



Chateau et al. 



(120081). 
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The main result of these studies is that, when model materials are carefully designed to corre- 
spond to the theoretical case of homogeneous and isotropic distribution of monodisperse hard 
spheres interacting only mechanically through a yield stress fluid, there is a theoretical relation- 
ship linking the dimensionless linear properties of such suspensions to their nonlinear properties 
that is in very good agreement with the experimental data. As regards their elastoplastic prop- 
erties, this relationship reads 



TM)/rM = V(l-'/')G'(0)/G"(O), (1) 

Tc being the yield stress, G' the elastic modulus, and (j) the particle volume fraction. They also 
observed that both the elastic modulus and yield stress are monotonically increasing functions 
of the solid volume fraction (j) which seem to diverge when (f) tends toward 0.57. Furthermore, 
their data are well fitted to a Krieger-Douglierty like equation 

G'(0) = G'(0) X (1 - (/-/^rn)"^-^"^™ with ^rn = 0.57 (2) 

and to its nonlinear generalization obtained by putting Eq. [2] into Eq. [1] 



tM) = rc(0) X ^(l-<^)(l-</,/^^)-2.5<^™ ^ith (f)^ = 0.57 (3) 

In this work, we again restrict to situations where there is a scale separation between the paste 
microstructure and the noncoUoidal particles in suspension and we focus on the purely mechanical 
contribution of the parti cles to the pas t e beha vior. For this purpose, we use the experimental 



procedures described in 



Mahaut et al 



(|2008al ). Accordingly, the experimental investigations 
reported in this paper focus on the behavior of the pastes in their solid regime, i.e., on the 
influence of particles on the elastic modulus and yield stress. 

In Secjlll we briefly present the materials employed and the experimental setup. Elastic 
modulus and yield stress measurements are shown in Sec. IIIII Then, in Sec. IIVI we present a 
packing model and show that it allows to account for the dependence of the rheological properties 
on the particle mixture composition. Finally, in Sec|Vl we combine the packing model estimates 
with the experimental data, and we demonstrate that it is possible to accurately predict the 
overall properties of the studied suspensions. 
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Figure 1: Shear stress vs strain when slowly shearing the material from rest at 10~^ for 
a pure emulsion (squares) and for a bidisperse suspension of 30% of 80 /xm polystyrene beads 
mixed with 70% of 315 /um polystyrene beads with cj) = 0.30 in the same emulsion (stars). 



II Materials and methods 



In this section, we briefly present the main features of the suspensions we studied and of the 
experimental procedures we used for the elastic modulus and yield stress measurements. It is 
recalled that these procedures were designed to study the purely mechanical contribution of an 
isotropic di stribution of particles to the yield stress fluid behavior. The interested reader is 



referred to 



Mahaut et al 



(|2008al )] for a more detailed presentation. 



A Materials 

In order to evaluate the purely mechanical influence of the particles on the behavior of the 
paste, we designed materials to ensure scale separation between the matrix (the yield stress 
fluid) and the partic 



Mahaut et al. 



es. We chose to use only inverse emulsions as a suspending fluid because 
(|2008al ) obtained the greater stability and the more reproducible results with this 
material in the monodisperse case. Moreover, the inverse emulsion behavior is very close to the 
ideal elastoplastic behavior of a yield stress fluid (see Fig. [T|). The emulsion is a water-in-oil 
emulsion which microstructure scale is given by the droplets size, of order 1 /xm from microscope 
observations. As the continuous phase, we use a dodecane oil in which a Span 80 
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emulsifier is dispersed at a 7% concentration. The dispersed phase is a 300 g/1 
CaC12 solution dispersed in the oil phase at 6000 rpm for 1 h. 

The particles are spherical beads. We used either polystyrene beads of density 1.05 g/cm^, or 
glass beads of density 2.5 g/cm'^. We used bidisperse mixtures of polystyrene beads of diameter 
80 nm and 315 /im (particle size ratio A = 3.94) and bidisperse mixtures of glass beads of 
diameter 40 and 330 fj,m (A = 8.25). The bimodal spherical beads are mixed together 
in the dry state. The particle mixture composition is defined by the fine particle proportion 
= Vp/iVp + Vc) where Vp (resp. Vc) denotes the fine (resp. coarse) particle volume fraction. 
The particles are then dispersed in a volume Vi of the suspending fiuid so that their total volume 
fraction is = (Vp + Vc)/{Vp + Vc + Vg). Finally, the fiuid-particle mixture is stirred manually in 
random directions in order to homogenize it and obtain an isotropic material. In all cases, the 
paste yield stress was sufficient to avoid particle sedimentation in the solid regime. The critical 
conditions for which a spherical object would fall under the action of gravity through 
a yield stress fluid at rest is obtained from the balance between the gravity force, 
the buoyancy and the drag force. A sphere of radius a immersed in a yield stress 
fluid with yield stress Tc will not move under the action of gravity if 4/3(/9 — pc)gT^a? 
is lower than ll k rira^Tr with kc ~ 3.5, pc the yield stress fluid density and g the 
gravity 



CoussotI (j2005l )]. In our experiments, the critical drag force is at least 30 
times greater than the gravity force, which ensures stability at rest (thus for the 
elastic modulus measurements). When the suspension is sheared at low shear rate 
, shear-induc e d sed imentation may occur with a velocity V 2/9(p — pcjga^ /{Tc/j) 



Ovarlez et al. 



ii]. o 



(|2010l )]. Our yield stress measurements (see below) are performed 



at very low shear rate (7 = O.Ols^): this yields V < 0.1/xm.s^^ for all the studied 
materials, i.e. shear-induced sedimentation can be considered as negligible over the 
duration of the experiments (less than 100 s). 

B Rheological methods 

The experiments were performed within a vane in cup geometry (inner diameter di = 25 mm, 
outer cylinder diameter de = 36 mm, height H = 4:5 mm) on a commercial rheometer (Bolilin 
C-VOR 2000). In order to prevent slippage at the walls, we used a six-blade vane as an inner 
tool immersed in a outer rough cylinder of roughness size equivalent to the size of the largest 
particles. We measured the elastic modulus G'{4') through oscillatory shear experiments at a 
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single frequency in the linear regime, and the yield stress Tc{(f)) through a single measurement 
at a small constant velocity (7 = O.Ols^^) on each sample (see an example in Fig. [T]). The 
(static) yield stress is defined as the shear stress plateau in a shear stress vs. shear 
strain plot (Fig. [1]). It is worth noting that v^e carefully checked that the same 
yield stress is measured whatever the (low) rotational velocity chosen to drive the 
inner tool is (i.e. 7 < O.Ols^^). These methods where shown to provide fair estimates of the 
relative evolution of the elastic and yield stress properties with the part icle volume fraction, a nd 



to ensure that the studied suspensions are isotropic and homogeneous 



Mahaut et al 



(j2008a|)]. 



Ill Results 

In this section, we present the elastic modulus and yield stress measurements performed on all the 
systems we studied, and compare the results to what was observed in monodisperse suspensions. 
Elastic modulus and yield stress measurements of various monodispe rse suspensions of beads 



embedded in yield stress fluids were obtained by 



Mahaut et al 



(2008a). 



A Elastic modulus measurement results 

The results of the elastic modulus measurements performed on all the suspensions we designed 
are summarized in Tab. [1] and depicted in Fig. [2] as a function of the solid volume fraction (p. 
Two size ratio (A = 3.94 and A = 8.25 ) and several fine particle proportion values ^ were tested. 

As classically observed for polydisperse viscous suspensions of different compositions, the 
experimental points do not fall onto a single curve. Nevertheless, we observe that the dimen- 
sionless elastic modulus is an increasing function of the solid volume fraction (/) when A and ^ 
are given. Moreover, all the data points for bidisperse suspensions fall below the law Eq. [2] valid 
for monodisperse suspensions. These results are perfectly consistent with the findings of the 
literature described in the introduction. It is also worth notin g that all the experimen t al ela stic 



modulus points fall above the Hashin-Shtrikman lower bound 



Hashin and Shtrikman 



19631)]: 



G'(,/.)/G'(O)>(2 + 30)/(2-20) (4) 

which is a theoretical lower bound computed in the general case of a biphasic material (an 
infinitely rigid phase embedded in a linear elastic phase) isotropic both at the microscopic and 
the macroscopic scales. 
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Figure 2: Dimensionless elastic modulus G'{(f),X,(,)/G'{0) as a function of the particle volume 
fraction (p for all the studied materials (all particle size ratio A and fine particle proportion 
^). The solid line is the law Eq. [2l The dotted line is the Hashin-Shtrikman lower bound Eq. HI 
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Table 1: Dimensionless elastic modulus G' {(j), \,^)/G' (0) as a function of the particle volume 
fraction for all the studied materials (all particle size ratio A and fine particle proportion 
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Table 2: Dimensionless yield stress rc((/>, A, i^)/tc(0) as a function of the particle volume fraction 
for all the studied materials (all particle size ratio A and fine particle proportion ^). 

B Yield stress measurement results 

We now present the results of the yield stress measurements performed on the different suspen- 
sions we studied. The experimental data are gathered in Tab. [2] and graphically represented in 
Fig. 3 as a function of the solid volume fraction 0. 

The dimensionless yield stress exhibits the same trends as the dimensionless elastic modulus. 
It is an increasing function of the solid volume when the particle size ratio A and the fine particle 
proportion ^ of the mixture are give n. Moreover, similarly to what was observed for monodisperse 



suspensions 



Mahaut et al 



(j2008al )]. for a given particle mixture, the yield stress increase with 



the particle volume fraction is much lower than the elastic modulus increase. 



IV Divergence of the rheological properties: The rigidity thres- 
hold of contact network 

Although the data shown in Sec. Illll are scattered, we observe that all the G' and Tc curves have 
basically the same shape and tend to diverge for a value (pm of the volume fraction that seems to 
depend both on the size ratio A and on the fine particle proportion ^. As a first step towards the 
modelling of the rheological properties of the suspensions, it thus seems important to be able to 
predict the value of (pm- This is the purpose of this section. 
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Figure 3: Dimensionless yield stress Tc{(j), A, (,) /tc(0) as a function of the particle volume fraction 
(j) for all the studied materials (all particle size ratio A and fine particle proportion ^). The 
sol id line represents the law Tc(0) = Tc{0)\/ (1 — 0)(1 — 4>/(j)m)~'^'^'^"^ with (prn = 0.57 proposed 
by Chateau et al. ( 20081 ) for the monodisperse suspensions. 



A Packing model 



Mahaut et al. 



(l2008al ) have shown that the elastic modulus G' of an isotropic suspension of 
monodisperse spherical beads is well fitted to the Krieger-Dougherty like law Eq. [2] with 0^ = 
0.57. The quantity (pm for which G' diverges is most probably no more than the 
contact rigidity threshold of the suspension: if < (pm the particles interact (hy- 
drodynamically) through the suspending fluid. If > (prnj the particles are close 
enough so that there is contact network spanning the whole material and that in- 
terparticle interactions (contact, friction, elasticity of the beads material ...) play 
a major role, leading to the apparent divergence of G'. Interestingly, the overall yield 
stress of monodisperse suspensions also strongly increases when (p tends towards the same value 
(pm = 0.57. 

At this stage, 0.^ is just a fitting parameter. However, according to the observed trends 
of the overall properties of the bidisperse suspensions in Sec. IIIIIAI and Sec. IIIIBI there is 
strong evidence that this c ontact rigidity thre s hold is closely linked to the r a ndom close packing 



Chang and Powelll (119941 ): 



Gondret and Petit 



(Il997h : 



Probstein et al. 



Frankel and Acrivos 



19671 )]. This suggests to use packing models to predict how the value of the contact rigidity 
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threshold vary with the beads mixture composition. Fm'thermore, it is largely accepted that 
the packing density of an assembly of monodisperse beads depen ds upon the small gap inter- 
particle forces and the way the packing is produced. For example, iDong et al\ (j2006l ) obtained 
packing fractions ranging from to 0.64 depending on the conditions between particles (0.64 
for random close packing of frictionless cohesionless beads and for random loose packing of 
beads involving strong van der Waals forces) . When the cohesive forces between particles do not 
dominate over other forces, the random loose packing fraction , or more precisely, the sphere pack - 



ing at its contact r i gidity t 



Onoda and Liginer 



1990|); 



ireshold ranges fr om 0.5 to 0.58 



Dong et al. 



Cumberland and CrawfordI (119871 ) 



(|2006l )]. Then, the contact rigidity threshold (pm = 0.57 
determined experimentally for monodisperse suspension, assuming it is representative of our 
mixture process, may be used as a fixed parameter of a packing model. 

Various models have been proposed in the literature to address granular mixture problems. 
For a bidisperse beads mix, the packing density is a function of the particle size ratio A, th e 



fine pa r ticle 



Doddd (|l980h : 



p roportion £ a nd th e way the packing was obtained 



Ben Aim and Le Goffl (jlMZl) 



de LarrardI (jl999l )]. In most of the models of the literature, the packing process 
was taken into account through a scalar index. In this work, we do not need to model the 
mixing process because we always used the same procedure to prepare the materials. 
We will first compute the value of the contact rigidity threshold of our bidisperse 
packings taking the monodisperse contact rigidity threshold 0m as a free parameter; 



the value of 0^ will then be fixed in the sequel as the ex perimenta 
value 0.57 for our mixing process. We chose to use the model of 



de Larrard 



l y ob served 



(Il999h which 



aims at taking the geometrical interactions between particles of different size into account. In 
this model, two different configurations of bidisperse mixtures are distinguished: dominant fine 
particle configuration when coarse particles are embedded in a fine particle matrix and dominant 
coarse particle configuration when fine particles fill the empty spaces between coarse particles. 
Furthermore, two geometrical interactions are taken into account: 

Loosening effect : When one fine particle is inserted into a packing of coarse particles (which 
are thus dominant class) and if the fine particle is not small enough to locate in the empty 
space between the coarse particles, there is a loosening of the coarse particle packing which 
induces a decrease of the overall coarse particle density. 



Wall effect : When one large particle is inserted into a local packing of fine particles (dominant 



13 



fine particle configuration), the fine particle density is decreased near the coarse particle. 
This wall effect induces a decrease of the overall fine particle density. 



Bournonville et al 



(j2004l ) used a slightly modified version of the model of ide LarrardI (jl999l ) 



to estimate the packing density of dry polydisperse beads mixtures. They showed that it is 
possible to directly use the experimentally measured packing density associated with the par- 
ticular packing process used as a free parameter of the model. Moreover, they experimentally 
determined the function describing the loosening and wall effects described above. While they 
measured random close packing for monodisperse beads ranging from 0.59 to 0.62 for parti- 
cles with diameters between 50^m and 320/im, we used exactly the same equations to estimate 
the divergence threshold of bidisperse suspensions, the only difference being the monodisperse 
threshold value. Since we are interested only in bidisperse mixtures, we only give the relation- 
ships allowing to est imate the contact rigidity threshold for such a suspension. The interested 



de LarrardI (|l999l ): 



Bournonville et al. 



(2004)1 for a more detailed review of 



reader is referred to 
the model. 

To compute the contact rigidity threshold of a bidisperse mixture of particles of respective 
diameters dc and dp with dc > dp, one may first calculate the fine dominant particle density 
and the coarse dominant particle density (/>^*^°™ defined by 



iF-dom 



xc-dom 



i-(i-0 (i-'/'j; + 6fc(c-|| 



i-e(i 



(5) 



(6) 



where (f)^ 



(resp. </>^) denote the contact rigidity threshold of the monodisperse fine (resp. 
coarse) suspension, a^p (resp. hpc) a function describing the loosening (resp. wall) effect and 
^ is the fine particle proport ion defined above. In the sequel, we adopt <j)'^ = (t>Z, = 0.57 as 
determined experii nentally by 



experimentally by 



Mahaut et al 



Bournonville et al 



(l2008ar ). The functions Qcp and bpc were determined 



20041 ) . They only depend upon the particle size ratio 
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Figure 4: Contact Rigidity threshold 4>m of bidisperse suspensions as a function of the fine 
particle proportion ^, for several size ratio A = 2, A = 3.94, A = 8.25 and A = 20. 



A = dc/ds. 



1.13' 



0.57 



(7) 



1.79 



0.82 



(8) 



Finally, the bidisperse suspension contact rigidity threshold is defined by 



min(C"°°\C'°"^) 



(9) 



Model predictions are shown in Fig. [U where the contact rigidity threshold is plotted as 
a function of the fine particle proportion for several values of the size ratio A. 

Note that whatever the value of the size ratio, the maximum of the contact 
rigidity threshold is a singular point of the contact rigidity threshold versus fine 
particle proportion curve. Even if such a singular point v^as not experimentally 
observed for bidisperse packing of dry part icles, it seems to be the hallmark of 



theore ti cal model di s tingu ishing tv^o regimes 



(119991)]. 



Gondret and PetitI dlMZl); 



de Larrard 



de LarrardI (jl999l ) corrected this undesirable effect by introducing a scalar 



index accounting for the packing process. Other models have been built to predict 



15 





1 

■ 


i 1 i 1 i 


1 


1 ' 


1 


100- 


• 


>.=3,94;£=0,1 




















80 - 


▼ 


>.=8,25; £=0,1 










♦ 


X=8,25; £=0,2 






■ 






X=8,25; £=0,3 






f 


60 - 




A=5,25; £=0,5 












>.=3,94; £<>0,0; 0,1; 0,3 






» 


40- 






















20- 




± > ■ ' — 




J ; 


0- 


•- 










1 

0,0 


1 1 1 1 1 

0,2 0,4 


1 

0,6 


0,8 


1 

1,0 



Reduced volume fraction ((()/<t)^) 



Figure 5: Dimensionless elastic modulus G'{(j),X,(,)/G'{0) as a function of the reduced solid 
volume fraction (j)/(j)m for all bidisperse suspensions (all particle size ratio A and fine particle 
proportion The contact rigidity thresholds (j)m were calculated using the modified 
de Larrard model Eqs. [5] to [9l 



smooth variations o f the p acking density w^ith respect to the fine particle proportion 



Gondret and PetitI (jl997l )]. To our knov^ledge, it is not clear that smooth models 
are more accurate than non smooth ones. Furthermore, it w^ill be shown in the 
sequel that Eqs. [5] to [9] are accurate enough to satisfactorily predict the influence of 
polydispersity onto both the elastic modulus and the yield stress of our suspensions. 
Then, the existence of this singular point seems to be a detail at this stage, and 
does not pose any problem in the framework of this study. 

B Rheological properties vs Reduced solid volume fraction 

In Sec. IIVIIA.I we have modelled the value of the contact rigidity threshold (p^^ i-e- the volume 
fraction for which the suspension rheological properties should diverge. We now propose to plot 
the dimensionless elastic moduli G'{(j), X,£,)/G' (0) and the dimensionless yield stress Tc{(j)) / Tc{0) 
of all our bidisperse suspensions as a function of the predicted reduced solid volume fraction 
(^/(pm- These data are presented in Fig. [5] and [6j 

We first observe that, when plotted versus the dimensionless elastic modulus now falls 

onto a single master curve. As the experimental data for monodisperse suspensions also 
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Figure 6: Dimensionless yield stress Tc{(p, X,(,)/tc{0) vs the reduced beads volume fraction (j)/(j)m 
for all bidisperse suspensions (all particle size ratio A and fine particle proportion ^). 
The contact rigidity thresholds (f)^^ were calculated using the modified de Larrard 
model Eqs. [5] to O 



fall onto this curve, this curve should be fitted by a Krieger Dougherty like equation. 
This question is addressed in the sequel. It is not obvious that such a property holds for 
the dimensionless yield stress, even if the experimental data are less scaterred than when they 
are plotted as a function of the solid volume fraction. We observe that both the rheological 
properties tend to diverge for (j)/(j)m = 1- This last observation is a strong indication that the 
packing model we use is able to predict correctly the impact of the particle size distribution on 
the value of the volume fraction for which the rheological properties diverge. 



V Analysis and discussion 

In Sec. IIVI we have shown that plotting the dimensionless rheological properties of bidisperse 
suspensions as a function of the reduced solid volume fraction </>/0m allows to obtain curves 
diverging when the reduced solid volume fraction tends toward 1. The purpose of the present 
section is to provide close-form estimates for the rheological properties of the bidisperse suspen- 
sions. 
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A Elastic modulus 



We have observed that dimensionless elastic moduh of all bidisperse suspensions plotted as a 
function of the reduced solid volume fraction (/>/0m fall onto a single master curve. Then, by 
slightly modifying Eq. [2]so that the exponent does not depend on we obtain a new estimate 
for the elastic modulus of bidisperse suspension 



G'{(t>A.€)/G'ifi) = {l 



-1.43 



(10) 



which is shown as a solid line in Fig. [5j Satisfactorily, experimental data are very well fitted to 
Eq. [ini As we have 1.43 = 2.5 x 0.57 for the exponent value in Eq. I10|, the Krieger- 
Dougherty like Eq JlOl agrees with Eq. [2l However, even if a good agreement between the 
experimental data and the theoretical equation is obtained, it has to be noted th at first order 
series expansion of Eq. [10] with respect to the variable (j) yields the lEinsteinl (|l906l ) relation 



G'((/>A, 0/G'(0) = 1 + 2.50 + O2(0) 



(11) 



only for monodisperse suspensions, i.e. when (pm = 0.57. Discrepancy between estimates 
of t he viscosity of coii c entrat e d suspension of solid parti cles and Einstein' s equation is class! 



cal 



Chang and Powelll (Il9941 1: 



I 



Frankel and AcrivosI (119671 )]. Eilers formula 



Stickel and Powell 



()2005l )] complies with both the Einstein equation and high concentration limit, but of course. 



this equation does not write as a single function of the normalized solid volume fraction 
because Eq. [11] does not. As our experimental data are well fitted to the closed form estimate 
Eq. [To] we leave this problem aside in this work. 

In Fig. 7, we plot the experimentally measured dimensionless elastic moduli G'{4>X,^)/G'{0) 
as a function of the fine grain proportion ^, for fixed size ratio A = 3.94, for two solid volume 
fraction (j) = 0.4 and (p = 0.5. These experimental data are gathered in Tab. [3] 

Experimental G'((/>A, ^)/G'(0) shows a minimum at a value of ^ close to 0.40. This minimum 
should thus correspond to the optimal mixture of particles of size ratio A = 3.94. The experimen- 
tal data are also compared with Eq. [TO]in Fig. [7] the contact rigidity thresholds being computed 
thanks to the model presented in Sec. IIVI 

The agreement between the experimental data and Eq. [10] is rather good meaning that the 
granular packing model combined with the elastic Krieger-Dougherty equation capture the essen- 
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Figure 7: Dimensionless elastic modulus G'{(p, A, /G'{0) vs fine particle proportion ^ for suspen- 
sions of 80 fim and 315 polystyrene beads with solid volume fraction (p = 0.4 (filled triangle) 
and (j) = 0.5 (filled square). The dot (resp. solid) line is Eq. [TOl computed for (j) = 0.4 (resp. 
<j) = 0.5). The contact rigidity thresholds were calculated using the modified de Larrard 
model Eqs. [5] to [H 



tial features of the studied system. However, the theoretical dimensionless elastic modulus has 
a minimum at ^ = 0.30 while the smallest elastic modulus was measured for suspension with ^ 
close to 0.4. This discrepancy possibly comes from the packing model which was validated on dry 
beads packing experiments, whereas beads are suspended in a fluid in this work. Finally, close 
inspection of the data suggest that the dimensionless elastic modulus vs fine particle proportion 
curve is smooth while the tangent to the theoretical curve is discontinuous at the point where 
G'{(f), X,^)/G' {0) reaches its minimal value. As pointed out in Sec. IIV|, this singularity in the 
theoretical curve clearly originates from the packing model. Even if the packing model could be 
improved to obtain smooth theoretical curves, it is believed that the model used in this work is 
accurate enough given the dispersion of experimental data. 

B Elastic modulus vs yield stress 

Proposing estimates for the overall prope rties of a suspension of particles dispersed in a non- 



Newtonian fluid is challenging. However, 



Chateau et al 



(120081 ) have recently shown that it is 



possible to relate the overall yield stress of a suspension of particles isotropically dispersed in a 
yield stress fluid to its overall elastic modulus in its solid regime, provided that the heterogeneities 
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Elastic Modulus 


(/> = 0.4 


(/> = 0.5 


C = 0.0 


5.60 


21.01 


^ = 0.1 


4.68 


15.16 


C = 0.2 


4.56 


10.90 


C = 0.3 


3.87 


9.05 


C = 0.4 


3.84 


8.07 


C = 0.5 


4.16 


9.23 


^ = 0.6 


4.51 


8.71 


C = 0.7 


4.25 


11.45 


C = 0.8 


5.03 


13.78 


^ = 0.9 


5.43 


16.71 


C = 1.0 


5.26 


22.91 



Yield stress 


(/) = 0.4 


= 0.5 


^ = 0.0 


1.91 


3.33 


^ = 0.1 


1.84 


2.93 


C = 0.2 


1.81 


2.61 


C = 0.3 


1.94 


2.33 


^ = 0.4 


1.61 


2.14 


C = 0.5 


1.79 


2.58 


^ = 0.6 


1.77 


2.62 


C = 0.7 


1.81 


2.69 


^ = 0.8 


1.80 


2.92 


^ = 0.9 


2.04 


2.99 


C = 1.0 


1.94 


3.26 



Table 3: Dimensionless elastic modulus G'{(j), A, ^) /G'(0) and dimensionless yield Tc{(p, A, ^)/tc{0) 
stress as a function of the fine particle proportion ^, for particle size ratio A = 3.94. 

of the strain rate field over the suspending fluid domain can be neglected (see Eq. [TJ. To obtain 
this relationship, it was assumed that the particles are rigid, noncoUoidal and that there are 
no physicochemical interactions between the particles and the pas te. It is reca ll ed tha t our 
experimental procedure was designed to fulfill all these hypotheses 



Mahaut et al. 



(12008aD]. In 



this case, the apparent viscosity r]'^PP of the suspension sheared at a macroscopic 
shear rate 7 reads 



(12) 



where 1]^^^ is the apparent viscosity of the suspending fluid sheared at an effective 
strain rate 7eff that accounts for the shear rate locally experienced by the pure 
suspending fluid, and G(0) is the elastic modulus of the same suspension of particles 
(i.e. same size and same position) dispersed in an elastic material. Eq. 1121 simply 
says that in a nonlinear (non-Newtonian) medium, the viscosity increase of the 
suspension with the particle volume fraction (p is simply the same as in a linear 
(Newtonian or elastic) medium whose viscosity would be the apparent viscosity 
Vo^^iieff) of the sheared suspending fluid; in this case, an additional dependence 
of the apparent viscosity on comes from the 7eff dependence on 0. As the solid 
particles do not deform, the shear correction factor 7eff/7 is greater than one. It can 
be shown in the framework of a rigorous u pscaling approach to this problem that 
the optimal estimate of this quantity reads 



Chateau et al. 



mm] 



7eff, 



■/7 = - </')G(0)/G(O) 



(13) 
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in which the quantity 1 — (j) appears because the shear rate experienced by the 
suspending fluid is linked to the overall shear rate by means of an average equation 
over the interstitial fluid only. This explains why nonlinear properties cannot depend 
simply on the reduced solid volume fraction <j)/(t)m- 

This approach allows in particular predicting the value rc((/>)/rc(0) of the dimen- 
sionless yield stress of suspensions of noncolloidal particles in yield stress fluids. 
The apparent viscosity of a perfect plastic yield fluid indeed reads rc/7. Putting this 
equation into Eq. 1121 combined with Eq. 1131 finally yields Eq. [H 

This analysis still allows to physically explain why the relative yield stress in- 
crease is lower than the relative elastic modulus increase. Indeed, while the relative 
apparent viscosity of the suspension is equal to its relative elastic modulus, the shear 
rate experienced by the suspending fluid is greater than the overall shear rate pre- 
scribed to the suspension. As the apparent viscosity of a perfect yield stress fluid is 
a decreasing function of the shear rate, the localization of the shear rate lowers the 
consolidating eff'ect of adding solid particles. 

We have plotted in Fig. [8] the dimensionless yield stress rc(0, A, C)/'Tc{^) — 1 as a function of 
the dimensionless quantity y^(1 — (f))G'{(p, A, ^)/G'(0) — 1 in logarithmic coordinates for all the 
systems we studied in order to check that Eq. [T]is still valid. We observe an excellent agreement 
between our experimental results and the micromechanical estimate Eq. [1] (which is plotted as a 
straight line y = x in these coordinates). 

These results show that the experimental data are consistent with the hypothesis that a 
uniform estimate of the strain rate over the suspending fluid domain allows to accurately estimate 
the overall yield stress of a bidisperse suspension from its overall elasticity. 

C Yield stress 

Putting the estimate Eq. [10] for the overall linear properties of the suspension into Eq. [1] yields 
the estimate for the overall yield stress of the suspension: 

TcicP, A, O/Tc = V(l-0)(l-0/0„)-l-43 (14) 

We have plotted in Fig. [9] the experimental dimensionless yield stress Tc{(f), X,£,)/tc{0) — 1 as a 
function of the dimensionless quantity — — (j)/(j)m)~^'^^ — 1 in logarithmic coordinates 
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Figure 8: Dimensionless yield stress rc((/>, A, ^)/rc(0) as a function of y^(1 — (j))G'{(j), A, /G'(0) 
for all the bidisperse suspensions studied (all particle size ratio A and fine particle propor- 
tion ,^). The figure's coordinates were chosen so that the y = x line represents the theoretical 
relation Eq. [TJ 

for all the systems we studied in order to check that Eq. [Inaccurately describes the experimental 
data. 

Even if the fit is not perfect, it is believed that Eq. [T4]is accurate enough to predict the yield 
stress of the bidisperse suspensions studied in this paper. 

In order to compare the accuracy of Eq. [14] as an estimate of the overall yield stress of the 
suspensions with the accuracy of the elastic modulus estimate Eq. [TOl we have plotted in Fig. [TU] 
the experimental dimensionless elastic modulus as a function of (1 — 

At first sight, both estimates seem to accurately fit the experimental data. To quanti- 
tatively assess this point, we have computed the least square error for both the 
dimensionless yield stress (Fig. [9]) and dimensionless elastic modulus (Fig. [T0|) . For 
quantity y the least square error reads 



with N, the number of measured points, y^^^, the ith experimentally measured point 
and y^^ the associated theoretical value. We found Errc = 0.0071 for the elastic 
modulus estimate and Err^^^ = 0.018 for the yield stress estimate. Both computed 




(15) 



i=l 
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Log(((l-(t,)(l>/<t>J-''''r-l) 



Figure 9: Dimensionless yield stress Tc{(p, A, ^) /tc(0) as a function of — — (j)/(j)m)~^'^^ for 
all the bidisperse suspensions studied (all particle size ratio A and fine particle proportion 
^). The figure's coordinates were chosen so that the y = x line represents the theoretical relation 
Eq.m 




Figure 10: Dimensionless elastic modulus A, ^)/G'(0) as a function of (1 — (f)/ (j)m)~^'^^ for 

all the bidisperse suspensions studied (all particle size ratio A and fine particle proportion 

The figure's coordinates were chosen so that the y = x line represents the theoretical relation 
Eq.M 



23 



errors being of the same order of magnitude, it can be concluded that the yield 
stress estimate Eq J14l is quite as precise than the elastic modulus estimate Eq. IIOL 



VI Conclusion 



We have studied the elastic modulus and yield stress of an isotropic bidisperse suspension of 
noncolloidal particles in a yield stress fluid. We focused on the purely mechanical contribution 
of the noncolloidal particles to the ove rall properties o f the y ield stress fluid. To do this, we used 



materials and procedures designed by 



Mahaut et al. 



(|2008al ) to study the case of monodisperse 



suspensions. We observed that, as it is classically observed for particles suspended in a Newtonian 
fluid, the elastic modulus and yield stress of bidisperse suspensions are lower than the same 
quantities measured for monodisperse suspension of same s olid volume fra c tion. We showed 



Chateau et al. 



mm hnking the 



that the equation rc(0. A, O/rc(0) = ^{l - (/>)G'((/))/G'(0) of 
overall yield stress, the overall elastic modulus and the solid volume fraction still applies when 
bidisperse suspensions are considered. Additionally, we showed that the eflFect of the particle 
size heterogeneity onto the overall rheological properties of the yield stress suspension can be 
described by means of a packing model developed to estimate random loose packing of assemblies 
of dry particles. We have finally proposed closed form estimates for both the elastic modulus 
and the yield stress. This shows that it is sufficient to determine (pm and the dependence of the 
rheological properties on (j) to predict the behavior of bidisperse suspensions; this should remain 
true for more complex polydisperse cases. An extension of this study to the cases of anisotropic 



particle distribution and more complex polydisperse suspensions is planned for the future. 
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